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A 


STUDY  OF  REGIONAL  PHASES 


1. 


Propagation  and  excitation  of  regional  phases  in  western  Europe 
Recall  of  former  work. 


-■'''in  a  first  approach('i-98&) ,  we  have  evaluated  the  attenuation  due  to 
geometrical  spreading  and  anelasticity  of  recorded  regional  phases  for  a  set 
of  selected  earthquakes  which  occurec  in  France  or  around  the  French  terri¬ 
tory.  (Nicolas  and  al.  BSSA  vol.  72  n°6  pp.  2089-2106  December  1982)  "‘This 
purely  experimental  study  leaded  to  the  conclusion  that  a  global  attenuation 
factor  in  the  range  of  distances  from  100  to  1000km  could  be  defined  on  broad 

— *v 

band  data  (0.5  -  16Hz).  Global  attenuation  coefficients  v  (A  a  D  )  around  2 


for  ?_  and  S,.  ,  and  around  2.3  or  2.5  for  ?  and  L  were  obtained, 
n  n  g  g 

Assuming  an  attenuation  due  to  geometrical  spreading  o  around  0.5 
-o 

(A- a  (D  sm  D)  )  we  could  compute  the  anelastic  attenuation  coefficient  K 


(ae  ).  For  L  waves  propagating  in  western  Europe,  K  is  of  theorder. 
of  0.2  *  0.1  deg  ~ 1  to  compare  to  0.07  deg  *  for  Eastern  US  and  0.5  deg 


for  northern  USSR  ,  or  0.35  deg  for  the  south  of  Casman  sea  (Nuttli). 


Its  was  also  checked  that  the  quality  factor  Q  was  depending  on  fre¬ 
quency  for  fach  of  these  phases  Pn  ,  Sn  ,  P  ,  L  .  Even  if  an  .accurate 
computation  was  not  realistic,  it  was  possible  to  estimate  Q  around  some 
100  to  some  1000  when  the  frequency  increases  from  1Hz  to  16Hz  . 

—-“In  a  second  approach-f-1-962)  ,  we  tried  to  model  regional  phases  by 
using  the  Bouchon’s  method,  in  order  to  understand  the  influence  of  the 
source  (depth,  source  mechanism)  and  of  the  propagation  model  on  the  seis¬ 
mograms  built  at  various  distances.  The  Bouchon’s  method  has  the  advantage 
of  a  full  waves  method,  which  computes  exact  seismograms  (Campillo  and  al 
BSSA  vol.  74  n°l  pp.  79-90,  February  1984). 

In  addition  this  theorical  study  brings  an  estimation  of  the  geometrical 

attenuation  with  epicentral  distance  more  accurate  than  the  previous  one 

for  Lg  waves  of  the  form  p-0.83  .  for  waves  of  the  form  D  . 

The  excitation  of  P  wave  is  found  insensitive  to  the  deml.  of  the  source 

6 

within  the  crust  while  the  wave  amplitude  is  found  50  per  cent  higher 

for  a  source  in  the  upper  and  middle  crust  than  in  the  lower  crust.  The 


amplitudes  of  these  two  phases  drasti  callv  decrease  when  the  source  is 

below  the  Moho,  which  underlines  the  important  role  of  wave  guide  played 

by  the  crust  for  the  propagation  of  ?  and  L  .We  also  find  that  the 

E  •  § 

geometrical  attenuation  of  P  and  L  waves  is  indeoendent  of  source 

E  E 

mechanisms.  In  the  case  of  an  explosion ,  the  excitation  of  ?  is  insensitive 

g 

to  the  source  depth.  The  L  wave  amplitude  is  small  in  comparison  to  ? 

E  g 

and  L  waves  and  depends  on  the  clo^ness  of  the  source  to  an  interface 
E 

or  to  the  free  surface. 


Research  and  development,  during  19  84  by  data  analysis  of  Lg  phases 
Taking  into  account  the  last  results  obtained  on  the  attenuation  of  L 

C 

waves  due  to  geometrical  spreading,  we  have  made  an  attempt  to  improve 
our  former  evaluation  of  the  anelastic  attenuation  and  consequently  of  the 
Q  factor.  The  method  we  have  used  also  leads  to  the  evaluation  of  the 
source  spectrum  of  the  quake  and  of  the  station  response. 

The  L  phases  which  are  made  of  multireflected  S  waves  within  the 
E 

crust  (Bouchon  1982)  carried  on  informations  on  the  source,  the  propagation 
and  the  reception.  In  the  frequency  domain,  the  amplitudes  are  expressed 

as  the  product  of  these  three  terms  :  Source,  Attenuation  and  Reception. 

We  have  separated  the  influence  of  these  three  terms  by  studying 
spectra  recorded  on  the  french  seismic  network  of  LDG  . 

This  network  of  large  aperture  (800km  x  800km)  is  well  adapted  to 
that  study.  The  telemetered  data  from  its  28  SP  seismographs  (vertical  com- 
ponant)  are  sampled  at  50s.  per  sec.  and  recorded  on  digital  form.  Consequent 
ly,  spectra  study  from  0.5Hz  to  15Hz  is  possible  with  a  reasonable  definition 


:)  Data_anci__Data  processing 


L  waves  records  of  IS  earthcuakes  located  around  France  (Fig.l)  with 
E 

magnitudes  ranging  from  Mi^  =  3.2  to  =4.8  are  selected  (Table  I). 

Propagation  paths  range  from  150km  to  1200km  but  for  one  quake  the  maximum 

recording  distance  is  generally  two  or  three  times  the  minimum  distance  due 

to  the  dynamic  range  (feOob)  of  the  digital  recording.  Each  L  waves  tram 

E 

is  divided  in  three  parts  according  to  their  group  velocities  (Fig.  2). 
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6. 


Lgl 

=  3.6 

km  /s 

>  V 

\ 

/ 

3.1 

km  Is 

maximum  amplitude 

T 

"g2 

=  3. 1 

km  Is 

>  V 

> 

2.6 

km  Is 

beginir.g  of  coda 

T 

~g3 

=  2.6 

km  Is 

>  V 

> 

2.3 

km  Is 

end  of  coda 

For  each  L  part,  the  frequency  spectrum  is  computed  (pig- 3)  within  0.5 
an  d  1 5H  z . 


Each  amplitude  .V  (f,d)  for  the  i  station  and  j  quake  is  expressed 
as  follows 


Aij(f.d)  =  S. (f )  *  EXG(d)  *  AI(f,d)  *  FTS.(f) 


M  \ 

k  1; 


wnere 


-  £.(f)  is  the  source  snectrum  of  the  auake  i 

1  *  ■*  J 


EXG(c)  is  the  attenuation  term  for  geometrical  spreading. 


-  Al(f.d)  is  the  term  for  anelastic  attenuation. 


FTS.  (f)  is  the  response  function  of  the  station  i  . 

:  is  frequency  d  is  distance 

Ijfd 

with  AI  ( f ,  d )  =  e  Q  V 

V  =  group  velocity 
Q  =  quality  factor 

The  geometrical  attenuation  term  has  been  found  previously  of  the  form 
EXG(d)  =  d~°'  6 

Consequently  equation (1)  is  : 

log  Aij(f.d)  -  0.8  log  d  =  log  S.(f)  -  -i4£-  -  log  FTS.(f)  (2) 

1  W  v  1 

which  is  solved  by  iteration  process. 


a)  For  each  frequency  and  each  quake  STf)  and  Q(f)  are  estimated 
by  least  squares  (FTS.(f)  is  put  to  1). 

b)  For  the  whole  set  of  earthquakes  u <e  build  up  a  relation 

Q  =  Q  f b 

o 

c)  For  each  station,  we  estimate  the  perturbation  to  bring  to  its 
transfer  function  as  the  mean  value  of  the  perturbations  given  by  ah  the 
earthquakes  in  this  station ,  with  the  additive  hypothesis 

Z  FTS.(f)  =  1 

i  =  1  1 

this  process  is  stable  and  convergent  within  few'  iterations. 


2°)  Main_  results 

a)  Transfer  function  at  the  iem  station. 

3y  taking  into  account  the  station  response,  the  Q  factor  estimation 
for  is  significantly  improved  and  the  data  dispersion  is  reduced.  The 

transfer  functions  for  the  23  studied  stations  correspond  essentially  to  atte¬ 
nuation  or  amplification  by  a  factor  up  to  3  for  high  frequencies  (f>5Hz)  Fig.  4  . 

These  transfer  functions  can  be  explained  by  important  variations  of 
Q  factors  associated  with  superficial  layers  beneath  the  station  area.  They 
are  well  correlated  with  variations  of  the  seismic  noise  spectrum  recorded  in 
the  station,  precisely  for  high  frequencies  which  seems  to  confirm  their  super¬ 
ficial  origme. 

b)  Quality  factor 

It  has  been  computed  independantly  at  each  frequency  for  each  earth- 
auake.  After  removing  the  station  transfer  function .  the  data  lead  to  an  esti¬ 
mation  of  the  Q  factor  with  an  accuracy  of  the  order  or  better  than  30£  • 

For  ull  the  quakes  and  the  three  groups,  the  Q  factor  varies 

significantly  with  frequency  following  a  relationship  as 


Q 


(Fig .  3) 


Figure  4  :  Transfer  function  of  each  station  (continue  line),  seismic  noise  sped  nun  (lolled  line) 


FREQUENCY 


11. 


we 

obtain 

for 

L 

gl 

Q 

=  290 

f0. 52 

for 

L 

S2 

Q 

=  295 

^0.55 

for 

L 

g3 

Q 

=  350 

,0.50 

These  similar  relationships  for  L  phases  and  their  coda  aeree  with  Aki 

g  - 

hypothesis  (  1980). 


c)  Source  spectrum 

Each  source  spectrum  is  obtained  as  the  mean  of  the  spectra  obtained 
in  all  the  stations  for  the  considered  earthquake,  corrected  from  each  station 
transfer  function,  and  from  the  propagation  or  attenuation  term.  The  anelastic 
attenuation  term  beeing  the  one  obtained  globally  as  described  in  b)  . 

The  spectra -"we  have  obtained  have  a  consistant  dispersion  probably  due  to 
a  too  global  propagation  term  which  does  not  take  into  account  differences 
between  the  propagation  structures. 

Nevertheless  they  own  the  main  general  features  of  the  Brune  source. 

The  L  phases  (L  )  and  their  coda  (L  and  L  )  lead  to  similar  source 
g  *  gl  S2  S3 

spectra  (Fig.  6).  In  the  high  frequencies  domain  of  the  spectra,  the  asymp¬ 
totic  decrease  well  defined  for  the  most  important  quakes  is  between 


■1.5 


and 


•2.25 


The  cut-off  frequencies  f  between  1.25Hz  and  10Hz  are  function 
of  the  magnitude  following  the  relationship 


log(f)  =  1.  98  -  0.  33  Ml 


On  the  other  hand,  the  low  freauencies  level  S  (0)  of  each  source 

o 

spectrum  is  a  function  of  the  magnitude  following: 

log  S  =  C  ♦  1.  5  M. 

true  for  L  ,  ,  L  „  ,and  L  waves. 

61  g2  83 

( with  CL  =  CL  -  0.38  ;  CL  =  CL  -  0.68) 

£2  gl  S3  gl 


Jog  amplitude  Jog  amplitude  *°g  amplitude 


log  frequency  (Hz) 


log  frequency  (Hz) 


Figure  6  L_  spectra  obtained  from  synthetic  seismograms 

...  m  >-  2  /  _  _  *■*  r _ _*  -n  _ _ i  ~  \ 


13. 


Ill 


3C)  Conclusion 


Spectra  of  phases  and  of  their  coda  recorded  in  western  Europe 
for  18  earthquakes  have  allowed  us  to  evaluate  the  influence  of  the  three 
essential  effects  wnich  model  them  :  source,  anelastic  attenuation,  and 
station  effects. 

The  source  spectra  we  have  obtained  have  the  main  general  charac¬ 
teristics  of.  the  5rune  source. 


The  aualitv  factor  Q  similar  for  the  L  ohases  and  their  coda 

g 

is  depending  on  the  frequency  as  we  knew  already  : 


Q 


300  f 


0.  53 


It  is  a  mean  quality  factor  for  western  Europe. 

A  transfer  function  for  each  station  is  estimated.  It  is  well  corre¬ 
lated  to  the  seismic  noise  spectrum  variations  versus  a  mean  spectrum  over 
the  network. 

This  transfer  function  mostly  significant  in  the  high  frequencies  domain 
is  probably  representative  of  the  superficial  layers  attenuation  under  the 
station  area. 


Research  and  develonment  in  1984c  Numerical  modeling  of  L  nhases 

- - -  - * - g- - 

Taking  into  account  the  previous  results  both  on  geometrical  spreading 
and  anelastic  attenuations  of  waves  in  western  Europe  we  have  tried  to 

interprete  the  Q  factor  of  these  waves  as  an  effective  measurement  of  QcJ 
mean  value  of  the  quality  factor  of  these  waves  in  the  crust.  The  crustal 
model  that  we  use  was  inferred  from  long  range  refraction  experiments  in 
central  France.  It  consists  of  a  4  layers  crust  overlying  a  mantle  half-space. 

For  P  waves  we  assume  Q  =  2  Q  .  The  calculations  are  done  using 

p  s  c 

the  discrete  radial  wave  number  representation  of  the  seismic  field  (Bouchon 
1981).  The  calculations  are  done  for  epicentra!  distances  ranging  from  180km 
to  700km  and  the  results  are  shown  in  Fig.  7  where  they  are  compared  to  the 
case  of  a  perfectly  elastic  crust.  The  source  considered  is  a  point  of  vertical 
strike-slip  dislocation  at  a  depth  of  10km, 

The  seismograms  represent  the  vertical  component  of  the  crustal 
response  to  a  sudden  step  function  dislocation.  The  frequency  range  extends 


Figure  7  :  Influence  of  attenuation  on  I,  phases  (versus  distance) 


:rom 


affects  the  wave  forms  and 


0  to  5Hz.  The  attenuation  strongly 
the  amplitudes  and  results  in  an  enhancement  of  ?  waves  and  Ravleieh 
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The  corresponding  wave  spectra  are  presented  in  pig.  £  for  a  time 

window  corresponding  to  group  velocities  between  3.8  km/'s  and  2.5  km/s  . 

Attenuation  of  the  high  frequencies  with  increasing  distances  is  clearly  visibl 

At  180km  ,  the  effect  of  attenuation  is  small  and  the  L  wave  spectrum  is 

g 

almost  flat.  This  differs  from  a  result  we  reported  earlier  (Campillo  and  al. 
1984)  and  which  was  due  to  an  incomplete  evaluation  of  the  wave-number 
series  for  large  horizontal  wave  numbers  at  high  frequency. 


The  numerical  results  obtained  are  compared  with  the  data  in  Fig.  9 

after  having  normalized  the  theoricai  spectra  so  that  the  mean  values  of 

observed  and  synthetic  spectra  at  180km  are  equal.  The  comparison  with 

the  data  shows  a  very  good  agreement  which  supports  the  interpretation 

of  L  waves  attenuation  as  a  shear-wave  crustal  attenuation, 
g 

We  have  also  looked  at  the  effect  of  source  depth  and  at  the  distri¬ 


bution  of  Q 


within  the  crust  on  the  attenuation  of 


L  waves  bv  rav  tracing 
g  ' 


All  the  subcritically  reflected  rays  have  been  drawn  (Fig.  10)  radiated 


by  sources  located  at  0  km  ,  10  km  and  29  km  depth. 


Clearly  the  L  waves  radiated  by  the  three  sources  sample  with  different 
weight  different  regions  in  the  crust.  A  deeper  source  will  sample  the  crust 
more  uniformly  than  a  source  located  in  the  upper  crust. 


We  also  have  traced  all  the  subcritically  reflected  rays  radiated  by 
a  10km  deep  source  and  leaving  the  source  in  three  different  ranges  of 
take-off  angles  :  0  -  60°  ;  60°  -  120c  and  120°  -  180°.  The  near  horizontal 

rays  do  not  sample  the  lower  crust  while  the  more  oblique  one  which  are 
associated  with  longer  travel  time  sample  the  crust  quite  uniformly  (Tig.  11). 


Consequently,  we  have  studied  the  sensitivity  of  waves  to  a  zone 

of  strong  attenuation  in  different  regions  of  the  crust.  A  zone  of  low  Q  in 

the  sedimentary  layers  or  at  shallow  depths  has  little  effect  on  the  wave 

train  with  respect  to  the  results  of  the  perfectly  elastic  case  (Fig.  12a  ). 

A  zone  of  low  Q  at  the  base  of  the  crust  on  the  contrary  gives  a  sharp 

attenuation  of  the  L  wave  coda  (Fig.  12b  )  which  was  insensitive  to  the 

g 

presence  of  highly  attenuation  sediments. 


Figure  8  :  Mean  source  spectrum  of  one  cuake. 


F i p ur e  10  :  L  ray  tracing  for  different  source  depths 
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Figure  12  :  Synthetic  seismograms  at  various  distances  from  the  earthquake. 
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Comine  back  to  the  three  different  time  windows  L  ,  L  ,  L  for 

gl  S2  ?  3 

L  waves  trains  recorded  in  France,  the  weak  increase  of  the  Q  factor 
S 

with  decreasing  group  velocity  indicates  that  Q  factor  of  the  lower  crust 
is  greater  than  the  mean  value  of  Q  in  the  crust.  Our  results  are  compa¬ 
tible  with  the  interpretation  of'  the  frequency  dependency  of  Q  as  a  tecto¬ 
nic  indicator  :  stable  regions  are  associated  with  Q  which  are  monotonic 
and  weakly  increasing  -with  frequency  w-hile,  as  pointed  out  by  Aki  (1980) 
active  tectonic  regions  show  a  minimum  of  Q  for  frequencies  around  0.5-lHz 
which  implies  a  stronger  increase  of  Q  with  frequency  as  predicted  by  the 
theoretical  works  of  Sato  (1984). 


IV  -  Future  work 

Excitation  and  propagation  of  regional  phases  and  particularly  of 
phases  have  been  studied  theonticalv  up  to  now  within  homogeneous  crustal 
layers . 

Similarities  between  synthetics  and  recorded  data  came  out  although 
complexity  of  most  of  the  phases  is  important. 

We  have  seen  that  by  using  a  more  complicated  propagation  model 
(several  layers,  one  sedimentary  layer)  we  could  approach  in  same  cases  a 
more  realistic  representation  of  the  different  phases. 

W'e  would  like  now  to  model  the  propagation  of  regional  phases  through 
heterogeneity  zones  such  as  large  fracture  zones,  non  planar  layers  etc.... 

This  should  be  done  by  an  extension  of  the  aiscret  waves  numbers 
method,  with  an  analysis  of  the  interface  effect  (Huygens  Principe). 
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STATION  CORRECTIONS  AND  COMMON  SIGNAL 


When  an  array  receives  a  seismic  signal,  produced  by  a  quake  some¬ 
where  beneath  the  array,  it  is  recorded  in  each  of  the  stations,  after  a  fil¬ 
tering  process  different  for  each  of  them  and  "polluted1'  by  signal  generated 
noise . 

The  aim  of  the  processing  we  will  describe  is  to  invert  this  set  of  signals. 
Under  ad-hoc  hypothesis,  it  is  possible  to  extract  from  the  records,  the  common 
signal,  the  station  corrections  and  the  noise  (mainly  the  coda  of  each  record). 

Data  used  in  the  analysis. 

As  a  data  processing  is  always  very  dependant  of  the  kind  of  data  usee, 
let  us  first  shortly  describe  the  records  which  will  be  processed  in  the  example 
given  below. 

The  signals  were  recorded  by  the  French  Seismic  Network  ;  28  stations, 

short' period  vertical  seismometers  (1Hz)  distributed  over  France  :  "Clusters" 
of  stations  (30  to  40km  apart)  the  clusters  themselves  are  at  a  distance,  one  to 
the  other.of  the  order  of  200  to  300km. 

The  events  used  in  this  analysis  have  to  be  relatively  far  away  from 
France  in  order  to  record  the  same  seismic  phase  in  all  stations.  In  the 
examples  below,  the  epicenters  are  at  a  distance  of  the  order  of  40°  from  the 
centre  of  the  network,  in  Central  Asia. 

We  believe,  the  distance  between  the  station  to  be  relatively  important 
in  this  analysis  :  not  too  small,  nor  too  large. 

Mathematical  Analysis 

At  first,  by  taking  into  account  the  time  of  arrival,  all  records  are 
phased.  The  station  corrections,  we  will  introduce,  deal  only  with,  the  ampli¬ 
tude  of  the  signal  and  not  with  a  possible  time-term. 
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We  will  suppose  tho^e  station  corrections  to  be  only  a  constant  for 
each  station  (kj  for  the  itr'  station).  This  could  be  considered  as  a  crude 
approximation  of  the  filtering  of  waves  by  the  geological  enviror.nement  of 
the  station.-  But  it  is  not  as  drastic  as  it  may  seem  :  it  is  oossible  to  work 
on  narrow  band  filtered  signals  and  to  apply  the  same  study  for  each  frecuen- 
cy  band. 


Let  us  call 


,th 


e<  .  the  n  sample  of  the  signal  recorded  in  the  i  "*■  station. 


Sj,  ,  the  ntn  sample  of  the  common  signal  we  are  looking  for. 

b;  ,  the  n^~  sample  of  the  noise  in  the  station  (mainly,  wave 

generated  noise  and  coda). 

We  have  to  solve  (in  k,s  and  b)  the  eauations 


ei,n  =  ki-£n 


D, 


i  ,n 


In  fact,  s  is  not  the  common  signal  which  we  suppose  the  existence  somewhere 
n 

beneath  the  array,  but  this  signal  filtered  by  the  mean  geological  structure  of 
the  array. 

Those  equations  are  strongly  under-determinated.  The  main  constrain, 
we  w'ill  add  to  solve  them,  is  that  everything  in  common  between  the  records 
has  to  be  put  in  the  common  signal.  It  means,  using  a  least  souare  approach 
that 

X  =  z  £Kn‘ 

1  n 

has  to  be  minimized. 

The  problem  may  now  be  solve  with  only  one  under-determination  left  it  is 
possible  to  multiply  every  k  by  the  same  factor  if  we  divide  the  common 
signal  s  by  the  same  constant.  This  fact  will  be  expressed  in  an  other  wav 
just  below. 

The  normal  equations 

J _ X 

3  kj 

_3 _ X 

5  sn 
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e. 

J  ,n 


s 

n 


n 


c 

"n 


v  > 
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-  1  i  i ,  n 
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the  last  of  them  express  that  the  common  signal  is  a  weighted  sum  of  the 

records  and  allows  us  to  eliminate  s  in  the  first  set. 

n 


k.  e.  e.  =  k.  I(  I  k.  e.  )  (I  k  e  ) 

1  J.n  i.n  i  n  .  3  j.n  \  £  £>n 


Introducing  the  cross  correlation  C  i , j  <■'£  the  signals  recorded  in  the  i1^  and 

.th  „  .. 
j  stations 


C  i,j 


n 


e. 

i,n 


e. 

J.n 


we  get 


(  Ik.  )  (  IC  i,j  k.)  =  k.  I  I  k.  k  Cj ,  £ 
j  j  j  3  1  J  £  3  £ 


or,  using  the  "bracket"  notation 


Clk  > 


<  k  i  C  Ik  >  , 

<  k,k  >  * 


Obviously  k  >.  is  one  of  the  eigenvectors  of  the  matrix  C  .  Its  module  is  not 
defined,  expressing  in  an  other  way  the  last  under-determination  we  pointed 
above . 

The  expression  X  ,  to  be  minimized,  can  be  written  as 


X  =  1  C. 


<  ki  Clk  > 

i ,i  irnr> 

i 

the  eigenvalue  associated  to  the  eigenvector  \  k> 


has  to  be  the  largest  one. 


This  analysis  is  very  similar  to  the  principal  component  analysis  in 


statistic.  It  could  be  possible  to  expend  the  all  set  of  records  on  the  basis 
of  all  signal-vectors  associated  to  each  eigenvectors  of  the  matrix  C  .  But 
we  only  see  a  possible  physical  explanation  for  the  first  one,  the  common  signal. 


Examples  : 

We  do  not  wish,  here,  to  present  neither  the  study  of  station  correc¬ 
tions  in  our  array  versus  azimuth  or  some  other  parameter,  nor  the  study  of 
a  seismic  swarm. 

We  will  just  give  as  an  example  of  the  processing  on  two  quakes  in  Asia. 
lc  :  1982  Dec.  5  50°K  B0°E  m  =  6. 1 

2°  :  1983  Oct.  fc  49. 9N  78.  8E  m  =  6.  0 


Those  two  events  were  chosen,  among  others,  for  their  large  magnitude  and  to 
be  rather  close  one  from  the  other.  The  focal  mechanisms  of  both  events  are 
very  similar. 

The  best  is  to  first  show  some  plotted  signals. 

The  figure  captions  are  : 


Fig.  1  Bottom  part  :  5  records  of  the  first  quake. 

The  total  duration  of  the  signals  represented  is  20  seconds.  The 
amplitudes  are  normalized. 

Upper  part  :  records  in  the  same  stations  for  the  second  quake. 


Fig.  2  :  Common  signals  for  both  events. 


Fig.  3  :  For  comparison  sum  of  the  phased  records  for  each  event.  The 

common  signal  is  a  weighted  sum  of  records.  The  signals  plotted  or. 
this  figure  are  sum  without  weighting  factors. 


Fig.  4  ;  Signals  and  Coda  :  for  both  events,  original  signals  in  two  stations 

and  the  computed  coda  :  signal  minus  common  signal  (designated 

by  C  at  the  beginning). 
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Discussion 


What  could  be  considered  as  an  advantage  of  a  least  square  procedure  is, 
that  it  always  gives  a  solution.  It  is  nevertheless  very  important  to  study  the 
quality  of  the  results. 

At  the  first  look,  the  differences  between  the  plain  sum  of  all  signals 
and  the  commun  signal  are  very  small.  Only  the  first  onset  of  the  records  are 
not  the  same.  But  this  is  rather  important.  A  large  part  of  the  information  on 
the  source  mechanism  is  contained  m  this  early  part  of  the  signal  (in  that  case, 
we  know  the  same  mechanism  to  be  very  short).  And  so,  our  process  does 
really  increase  the  ability  to  study  this  mecanism. 

We  also  computed  the  station  corrections  kj  for  various  time  window 

length  (the  results  plotted  above  are  for  the  largest  of  the  windows  used). 

As  long  as  the  first  two  seconds  of  each  record  are  in  the  window  the  k. 

1 

are  stable  (within  5%).  For  both  quakes,  the  differences  are  less  than  10%. 


Future  Work 

Obviously  we  are  only  at  the  beginning  of  our  work.  We  have  to  study 
the  results  of  our  data  processing  system  : 

On  a  seismic  swarm,  to  obtain  better  source  function  and  to  compare 
the  source  from  one  event  to  the  others. 

-  On  the  array  :  to  obtain  station  corrections  versus  azimuth,  to 
apply  the  corrections  obtained  on  a  large  event^on  the  records  of  smaller  ones.. 

-  On  the  study  of  the  coda  others  preliminary  results  let  us  think 

that  within  the  cauda,  it  is  possible  to  define  secondary  arrivals.  We  wish  to 
localize  their  origin. 


